Strong and complex electron-lattice correlation in optimally doped E^S^CaCi^Os+tf 
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We discuss the nature of electron-lattice interaction in optimally doped Bi2Sr2CaCu2 0s+(5 sam- 
ples, using isotope effect (IE) in angle resolved photoemission spectroscopy (ARPES) data. The 
IE in the ARPES line width and the IE in the ARPES dispersion are both quite large, implying a 
strong electron-lattice correlation. The strength of the electron-lattice interaction is "intermediate," 
i.e. stronger than the Migdal-Eliashberg regime but weaker than the small polaron regime, requiring 
a more general picture of the ARPES "kink" than the commonly used Migdal-Eliashberg picture. 
The two IEs also imply a complex interaction, due to their strong momentum dependence and 
their differing sign behaviors. In sum, we propose an intermediate-strength coupling of electrons to 
localized lattice vibrations via charge density fluctuations. 

PACS numbers: 74.25.Jb, 74.72.-h, 79.60.-i, 71.38.-k 



In the last few years, several angle resolved photoe- 
mission spectroscopy (ARPES) reports suggested a sig- 
nificant interaction of dynamic lattice distortions with 
electrons in the cuprate superconductors p], but 
the basic nature of the electron-lattice interaction is still 
unclear. 

On one hand, the Migdal-Eliashberg theory, a standard 
model of the electron-phonon interaction in solids, seems 
to provide a basic framework for explaining [5|, [6[ a key 
feature of ARPES data, i.e. the "kink" in the ARPES dis- 
persion 

UMM- In this theory, the ARPES kink appears 
at phonon energy, uo p . For excitation energy (u) smaller 
than long-lived quasi-particles (electron dressed by 
virtual phonons) form. For uu > uo pl short-lived electrons 
scatter strongly with real phonons. For uo ^> u p , the ef- 
fect of the scattering becomes negligible. This familiar 
Migdal-Eliashberg picture has been also generalized to 
cases where other bosonic excitations are involved @, Q . 

On the other hand, a large number of experimental 
reports jjj |, EE El El El, EI El, El, EE E3 suggest an 
interaction strength beyond the Migal-Eliashberg theory. 
Such a proposal is plausible also from the ARPES point 
of view, because as momentum approaches the Brillouin 
zone boundary, uo p 70 meV) becomes comparable to 
the effective band width governed by Fermi velocity [19[ , 
weakening the Migal-Eliashberg theory, and the actual 
quasi-particle weight Z measured by ARPES is rather 
small, - 0.1 [20[. 

To date, there has been no ARPES study reconcil- 
ing these two important points of view. This absence 
appears to be due to the complexity of the cuprate 
physics, which makes it difficult to disentangle the effect 
of electron-lattice interaction from the effect of strong 
electron-electron interaction. As we demonstrated re- 
cently [3], ARPES on isotope substituted samples is a 



unique method that can directly sort out the role of the 
lattice to shed light on this complexity. 

In this Letter, we provide evidence of an anomalously 
large isotope effect (IE) in the ARPES peak width for 
optimally doped Bi2Sr2CaCu208+5 superconductor, in 
addition to the IE in the ARPES peak position reported 
earlier [3]. Detailed comparison of the data with vari- 
ous theories of the electron-lattice interaction provides 
compelling evidence that the strength of the interaction 
is intermediate, i.e. beyond the Migal-Eliashberg regime 
but not in the small polaron regime. This and another 
feature, the local nature of the lattice vibrations involving 
only a few lattice sites @, [H| , are an important refine- 
ment of an anisotropic electron-lattice interaction model 
proposed before [ESQ- Finally, we believe that our 
results give important clues for understanding the nodal- 
anti- nodal dichotomy observed in the cuprates H, E2I as 
well as the interplay between the electron-electron inter- 
action and the electron-lattice interaction, as proposed 
recently jjj HE HI • 

We discuss data taken at low temperature (25 K), 
where the electron-lattice interaction is strongly en- 
hanced [3[ well below T* (pseudo-gap temperature) « T c 
(92 K for 16 O and 91 K for 18 O). Details of the experi- 
ment are described elsewhere [3]. Throughout this Let- 
ter, we use the term "high energy" or "low energy" to 
refer to the magnitude \E — Ep\ [Ep is the Fermi en- 
ergy), high or low meaning relative to the kink energy in 
particular. 

In Figure 1 we show raw ARPES data along two im- 
portant momentum space cuts, cut a (nodal) and cut b 
(off-nodal or near-anti- nodal [25[). A more detailed mo- 
mentum dependence is described elsewhere [3|, |26| for the 
peak position and below for the width [25]. The energy 
distribution curves (EDCs; intensity as function of energy 
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FIG. 1: Comparison of ARPES EDCs for three samples 
with different isotope treatments, normalized to same peak 
height; normal ( 16 0), isotope-substituted ( 18 0) and isotope 
re-substituted ( 16 Or) samples for cuts a and b, as indicated 
in the inset of panel b along with the Fermi surface. Data for 
each angle were integrated over a 1° angle window in favor 
of statistics. The inset saddling the two panels shows Migal- 
Eliashberg simulations for 16 O and 18 O, as described in text. 

at fixed momentum value) displayed here show clear IE 
for both the peak position and the peak width. Previ- 
ously [3], momentum distribution curves (MDCs; inten- 
sity as function of momentum at fixed energy value) were 
used to discuss IE in the MDC dispersion and a small 
red-shift of the kink energy. In this Letter, we go further 
to discuss the line width, which requires more care in 
analysis. For this reason, we discuss EDCs rather than 
MDCs and use the half width at half maximum on the 
right side (HWHM#) of EDC peak to quantify the peak 
width [27], as illustrated in panel a. As we will show, the 
EDC analysis gives the same result as the MDC analysis 
in terms of the peak position. 

Two observations can be added [3] to support the in- 
trinsic nature of the observed IE. First, high energy 
peaks in Figure 1 sharpen up upon isotope substitution, 
reinforcing the argument [3] that a simple band structure 
effect or a simple disorder effect, arising from a stronger 
disorder for the 18 O sample due to the substitution pro- 
cess, cannot be the origin of the observed IE. Second, the 
EDCs for the re-substituted sample (green) are nearly 
identical to those of the 16 O sample. 

It is quite clear that the Migal-Eliashberg theory (in- 
set) predicts a small IE overall, in strong contrast to the 
data. For the Migal-Eliashberg theory, we chose the fi- 
nite temperature formalism [5[ for electron coupled to 
seven Einstein phonon modes, equally spaced from 10 to 
70 meV, with the maximum determined by neutron scat- 
tering data (28[. For this model, we refer to this max- 
imum frequency as uo p . For the purpose of this paper, 
this choice of phonon frequencies is equivalent to a more 



realistic [28[ choice of a fewer number of high frequency 
modes, while doing a better job of describing a smooth 
kink [3[ seen in the data. We use a gap-less linear band 
with vp = 1-8 eVA for e(fc), and 0.5 for the dimension- 
less electron-phonon coupling constant A [5]. For 18 O, 
only cjp is modified, in accordance with the harmonic 
phonon model. We have checked extensively that use of 
additional features such as anisotropic coupling [29[, full 
band structure and superconducting gap [y, [7|, [8[ does 
not affect the discussion presented below. 

In Figure 2, we present a detailed analysis of the IE 
introduced in Figure 1. Panels a and b show "peak IE," 
namely the isotope dependence of EDC peak position. 
Likewise, panels c and d show "width IE," the isotope de- 
pendence of EDC peak width. Several important points 
can be made. First, both IEs are large at high energy, 
while small at low energy [30]. This result shows a com- 
plete agreement with the previous work 0] as the peak 
IE is concerned, and further shows that the kink energy 
(« the arrow position in panel b for both cuts) continues 
to separate small IE and large IE also for the width IE, as 
it separates the sharp coherent peak and the broad inco- 
herent peak [3[. Second, the Migal-Eliashberg simulation 
(light red lines in b,d) poorly explains the data in gen- 
eral and the data at high energy in particular. Third, the 
IEs are enhanced for the off- nodal cut, consistent with a 
strongly momentum dependent electron-lattice coupling 
[6[ and making the failure of the Migal-Eliashberg the- 
ory less severe for the nodal cut. Fourth, the sign of 
the width IE is unchanged going from cut a to cut b, 
i.e. the peak always sharpens up for 18 O, while the sign 
of the peak IE changes from cut a to cut b. These dif- 
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FIG. 2: (a,b) IE in the EDC peak positions for cuts a and 
b. (c,d) IE in the EDC widths (HWHM H ). The IE is defined 
as 16 O — 18 O, and note that the sign- reversed IE is plotted in 
panel b. The horizontal energy axis correponds to the peak 
position for the 16 O sample. Gray dashed lines are IEs for 
16 Or. For comparison Migal-Eliashberg simulation results 
are shown as light red lines in panels b,d. 
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FIG. 3: IE of the EDC width as a function of the super- 
conducting gap A, for six momentum cuts [25j] shown in the 
right panel. For each cut, low energy (-70 meV to meV) 
average and high energy (-250 meV to -70 meV) average of 
IE are shown with error bars. The line shown for high energy 
IE is guide to the eye. 



fering sign behaviors are difficult to understand within 
the usual self energy analysis scheme [31], because the 
Kramers-Kronig relationship connects the real part and 
the imaginary part of self energy [3l[ , determined mainly 
by the peak position and width respectively. This means 
that the self energy analysis scheme as it has been em- 
ployed to date under two key assumptions, momentum 
independent self energy and high energy data cutoff at 
~ 0.3 — 0.4 eV, is unreliable to the accuracy relevant here. 

In Figure 3, we summarize the momentum dependence 
of the width IE up to cut b [25]. The low energy IE 
is small, and vanishes within error bars. In contrast, 
the high energy IE is finite and shows an approximately 
linear correlation with the superconducting gap up to cut 
b [HI , as the peak IE Q . The corresponding IEs for the 
Migal-Eliashberg theory is negligible, 2 meV at the most. 

The failure of the Migal-Eliashberg theory can be 
traced back to the the single phonon loop approximation 
for electron self energy, subsequently resulting in a small 
IE (shown by light red lines in Figures 2b and 2d) gov- 
erned at the most by the isotope dependence of the single 
phonon energy. For instance, within this theory the high 
energy line width is given by ~ uo pi explaining the small 
IE at high energy (Figure 2d) as well as explaining why 
the ARPES width of this theory is by a factor of « 2 
too small in comparison to the observed width (Figure 
1). In the literature, there have been two different ap- 
proaches to remedy this failure of the Migal-Eliashberg 
theory. First, Seibold and Grilli [7( used a more general 
Migal-Eliashberg model which involves A that is strongly 
momentum and isotope dependent. The physical reason, 
especially for the latter, is that the boson that couples to 
electrons in this case is critical charge order fluctuations 
of correlation length of a few lattice constants, instead 
of phonons extended throughout the lattice. Using this 
generalization, it was shown that the model can qualita- 
tively explain the sign change seen in Figures 2a, b but 
not the IE seen in Figures 2c, d, predicting a broader peak 



for 18 O. Second, models [32|, [33[ incorporating interac- 
tion strength going beyond the Migal-Eliashberg theory 
were shown to be able to explain the large (small) iso- 
tope effect at high (low) energy and also the sharpening 
of ARPES peak induced by 18 O. On the other hand, the 
sign change of Figures 2a, b could not be addressed, pre- 
sumably due to the simple form of interaction assumed. 
It is significant to note that in these calculations the 
binding energy of peak increases while the peak width 
decreases, a non-trivial result in qualitative agreement 
with our nodal cut data (Figures 2a, c). 

Based on these results, we propose that a combination 
of these two approaches is required to be totally success- 
ful, as we further elaborate now. 

In Figures 4a, b, we show line shape simulations for cut 
b using the Holstein model in the strong-coupling (or the 
anti-adiabatic) limit where a simple perturbative solu- 
tion can be obtained [34]. For input parameters of the 
Holstein model, we use uo p = 70 meV (66 meV for 18 O) 
and tight binding fit of e(k) ^ to the 16 O ARPES data. 
The interaction parameter A was taken to be isotope de- 
pendent [7[, 0.9 for 16 O and 0.8 for 18 O, in order to give 
the correct sign for the peak shift. Note that the peak 
sharpening for 18 O remains true even if A is taken to be 
isotope independent. The essential characteristics of the 
line shape simulation is the strong multi-phonon "shake- 
up" features, which occur at harmonics of co p (panel a) 
and which we postulate will broaden, due to phonon con- 
tinuum and strong electron-electron interaction in solids, 
into a single peak (panel b), as was also suggested re- 
cently [13]. In this way, the simulation in panel b suc- 
cessfully reproduces the small IE for sharp low-energy co- 
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FIG. 4: (a,b) Simulation of EDCs at k — for cut b using 
a small polaron theory [34]. The line shape simulations in- 
clude a 10 meV (panel a) and an additional 100 meV (panel 
b) FWHM Lorentzian broadening, the latter only for multi- 
phonon peaks, (c) Proposed schematic diagram of ARPES 
kink for intermediate and weak electron-lattice interaction. 
The thickness of the line used for CP or IP dispersion crudely 
represents the ARPES line width. The intensity of CP or IP 
is not indicated, but may be roughly inferred from the inset, 
which corresponds to the main panel what panel a is to panel 
b. See text for more discussion. 
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herent peak (CP) and the large IE for broad high-energy 
incoherent peak (IP). As shown, IEs as large as ~ 30 
meV are generic features of this simulation. In addition, 
the magnitude of the ARPES line width is also realistic. 
These important results make an explicit demonstration 
of the multi-phonon physics proposed within a dynamic 
spin-Peierls picture [3[. 

The strong coupling theory used here, while capturing 
the essence of the multi-phonon (polaronic) physics, has 
two obvious shortcomings in explaining the ARPES kink, 
a virtually non-dispersive CP with a very small weight 
Z <C 0.1 and a completely non-dispersive IP [34]. Both 
these shortcomings disappear as the interaction strength 
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is reduced [36|, |37|, [38|, and thus it is suggestive that 



the interaction at optimal doping lie in the intermedi- 
ate regime, where there is still a significant multi-phonon 
contribution to the electron self energy and both CP and 
IP (especially the latter) show strong dispersions. While 
it may be naively expected that the IE will decrease as 
the interaction strength is reduced, studies [H, [IH show 
that actually at intermediate couplings the IE is anoma- 
lously enhanced, strengthening our qualitative argument 
here. 

Our results clearly call for a reconsideration of the 
usual ARPES kink picture based on the MDC analy- 
sis within the Migal-Eliashberg theory. In Figure 4c, 
we propose a schematic picture that is applicable to in- 
termediate coupling regime as well as Migal-Eliashberg 



The 



regime, based on numerical simulations [36 
inset shows numerical simulations (Figure 4 of Ref. 
the intermediate coupling regime, showing distinct multi- 
phonon branches reminiscent of multi-phonon peaks in 
panel a. Due to the additional broadening mechanisms 
discussed above, these multi-phonon branches will tend 
to merge into a single IP branch, as shown in the main 
panel. The dispersion of IP at high energy is similar, but 
not identical, to e(k) [H, [I?!. Within this picture, the 
kink anomaly still occurs at uo pi the onset of multi-phonon 
continuum (gray area), consistent with a small isotope- 
induced shift [3]. Note that this picture places EDCs 
as more basic quantities than MDCs. Subsequently, the 
kink energy is best defined as the separation energy be- 
tween CP and IP branches of EDCs, and the MDC dis- 
persion with a kinky crossover (dotted orange line) is 
merely a consequence of weight transfer between the two 
branches. 

To conclude, we have shown that IEs in ARPES peak 
width and peak position give important clues to the na- 
ture of electron-lattice interaction in optimally doped 
Bi2Sr2CaCu208+5- Theoretical [13] and experimental 
(2lj works further support the emerging picture in which 
a cooperation of local lattice distortions and charge/spin 
order fluctuations produces an interaction effect on elec- 
trons beyond the limit of the Migal-Eliashberg theory. 
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